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Abstract. - We report on a comprehensive study of the electronic and magnetic structure of 
the striped nickelate Las/sS^/sNiCk. The investigation is carried out using band structure calcu- 
lations based on density functional theory. A magnetic structure compatible with experiment is 
obtained from spin-polarized calculations within the generalized gradient approximation (GGA), 
whereas inclusion of a local Coulomb interaction in the LDA+C/ framework results in a different 
ground state. The influence of the various interaction parameters is discussed in detail. 



Stripe phases are currently the focus of an intense ac- 
tivity as they are observed in a wealth of systems. Most 
of them are layered oxides, including cuprates [1-4], nick- 
elates [5-9], and manganites [10-12]. In such doped Mott 
insulators the structure of these phases varies from system 
to system: the domain walls in which the holes are primar- 
ily located can run parallel to the lattice axes, or along the 
diagonal. Furthermore, the distance between them is ei- 
ther compatible with one doped hole per two domain walls, 
in the so-called half-filled stripes, or with one doped hole 
per one domain wall in the filled stripes, as discussed by 
Zaanen and Oles [13]. Those various structures are real- 
ized in different systems, e. g., half-filled vertical stripes 
in the lightly doped cuprate Cai.ssNao.i^CuCbCk, and 
diagonal filled stripes in I^-^Sr^NiC^. Stripe formation 
manifests itself in a wealth of experiments, in particular 
by the appearance of peaks in both unpolarized and po- 
larized neutron scattering experiments. Indeed charge and 
spin order in La2- a; Sr a ;Ni04 is characterized by the wave 
vectors Q c = 27r(e,e) and Q s = 7r(l±e, l±e), respectively, 
with e ~ x corresponding to a constant charge density of 
one hole/Ni ion along the diagonal stripe. Incidentally, the 
electronic states involved in the stabilization of the stripe 
phases form mid-gap bands; in the first case these bands 
are partially filled, sometimes even half-filled, while in the 
second case these bands are filled (for a recent discussion 
see Ref. [14]). From the point of view of describing these 
phases by means of microscopic models it is commonly ac- 
cepted that the one band Hubbard model is suitable for 




Fig. 1: (Color online) Magnetic structure of the Ni atomes in 
the NiO-planes of Las/3Sri/3Ni04. Bright (yellow) circles cor- 
respond to Nil/Ni2 sites, dark (blue) circles to Ni3/Ni4 sites. 
The representation refers to the experimental CDFAS phase, 
see the text for details. 



the cuprates, whereas the orbital degeneracy incorporated 
in a two-band Hubbard Hamiltonian plays a prominent 
role when studying nickelates. In both cases it turns out 
that the parameters of the model deeply influence the type 
of ground state (see, e. g., Ref. [15]), and it would be de- 
sirable to obtain an estimate of them from ab initio cal- 
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semi-core states 


valence states 


Sr 


4s, 4p 


5s, bp 


Ni 


3p 


3d, 4s, 4p 


La 


5s, bp 


4/, 5d, 6s, 6p 





2s 


2p, 3s 



Table 1: Choice of semi-core and valence states. 



culations. Indeed, when comparing a doubly degenerate 
Hubbard model to a more realistic Hamiltonian containing 
e g states [16], the former predicts bond centered diagonal 
half-filled stripes to be the ground state, whereas diago- 
nal filled stripes are promoted in the latter case. Besides, 
the relevance of the Jahn- Teller coupling was emphasized 
by Zaanen and Littlewood [17]. Therefore a band struc- 
ture calculation embracing all these aspects is expected to 
provide insight into the mechanism leading to the stripe 
formation. As compared to the cuprates, quantum fluctu- 
ations should play a smaller role in the nickelates where 
the Ni 2+ ions carry spin 1. The approximations involved in 
the band structure calculations hence are expected to have 
less impact. In this work we focus on the doping x = 1/3, 
where the wave vectors characterizing spin and charge 
modulations coincide as Q s , c = §(L1)- This case cor- 
responds to the most stable stripe phase with an ordering 
temperature of T ~ 240 K, as revealed by specific heat [18], 
transport [19], and optical conductivity data [20]. For a 
more comprehensive review of the experimental situation 
we refer to Ref. [16]. 

The crystal structure of the Laa-xSr^NiC^ compounds 
is based on the body-centered tetragonal K 2 NiF 4 structure 
with space group I4/mmm [21], as for the La2- a; Sr 2 ;Cu04 
system [22]. Therefore, the tetragonal unit cell comprises 
two La2- a; Sr a ;Ni04 formula units. The compounds essen- 
tially consist of NiC>6 octahedra forming layers parallel to 
the afr-plane, where neighbouring layers are separated by 
La/Sr ions. To be more specific, the magnetic Ni atoms 
are arranged on a square planar lattice and interact via 
the atoms located midway between them. By symme- 
try, there are two crystallographically inequivalent NiO 
planes, which, however, are closely related to each other. 
We have obtained the structural parameters for our band 
structure calculations by interpolation between the exper- 
imental values at compositions of x = and x — 0.6. For 
the composition of interest, x — 1/3, the lattice constants 
amount to a — 3.83 A and c = 12.68 A. While the La/Sr 
sites are located at (4e) Wyckoff positions with parameter 
z — 0.3617, the Ni and 0(1) sites occupy (2e) and (4c) 
Wyckoff positions, respectively. Finally, for the 0(2) sites 
we again have to deal with (4e) Wyckoff positions, where 
the parameter is z = 0.174. In order to account for the 
realized spin and charge modulations, we set up a super- 
cell of the tetragonal unit cell. Comprising altogether six 
La5/3Sri/3Ni04 formula units with 12 La/Sr sites, our su- 
percell is given by the lattice vectors (1, 1, 0), (3, 0, 0), and 
(0, 0, 1), with respect to the tetragonal lattice. In order to 
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Fig. 2: (Color online) Partial Ni 3d and O 2p DOS (per unit 
cell) for spin-degenerate LeL 5 / 3 STi/ 3 ~Ni04, as obtained by the 
LBA+U method. 

satisfy the La:Sr ratio of 5:1, we have to attribute 10 La 
and 2 Sr atoms to the 12 La/Sr sites, where the details of 
the distribution do not affect our further argumentation. 
As a consequence, we have four classes of inequivalent Ni 
sites, called Nil, Ni2, Ni3, and Ni4 in the following. Sites 
Nil and Ni3 are located in one of the inequivalent NiO 
planes, while sites Ni2 and Ni4, respectively, correspond 
to them in the other NiO plane. The positions of sites 
Nil and Ni3 are translated into that of Ni2 and Ni4, re- 
spectively, by the vector (a/2, a/2, c/2). The charge and 
spin pattern covered by these structural prerequisites is 
illustrated in fig.[TJ We have diagonal stripes of Nil (Ni2) 
sites, separated by domains of Ni3 (Ni4) sites, thus reflect- 
ing the coincidence of the spin and charge wave vectors. 

The electronic structure calculations presented in the 
following rely on density functional theory (DFT). Con- 
cerning the DFT implementation, we apply the Wien2k 
package, a state-of-the-art full-potential code with mixed 
lapw and apw+lo basis set [23]. This code has shown great 
capacity in dealing with the interplay of structural relax- 
ation, magnetic interaction and electronic correlations in 
very complex materials [24]. In each of our calculations, 
the charge density is represented by almost 23,000 plane 
waves and the mesh for the Brillouin zone integration com- 
prises 60 k-points in the irreducible wedge. Our choice 
of semi-core and valence states is summarized in table [TJ 
For the supercell setup described above, we address cal- 
culations within the pure generalized gradient approxima- 
tion (GGA) as well as LDA+J7 calculations to account for 
an additional local electron-electron interaction. We use 
the Perdew-Burke-Ernzerhof [25] parameterization of the 
exchange-correlation functional for the pure GGA, while 
the LDA+J7 calculations are based on the SIC scheme de- 
scribed in [26,27]. In each case, we use a band structure 
calculation where we artificially enforce spin-degeneracy as 
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Fig. 3: (Color online) Partial spin-majority and spin-minority 
Ni 3d DOS (per Ni atom) for spin-polarized La5/3Sr!/ 3 Ni04 in 
the ADFAS phase, as obtained by the LDA+f/ method. 



Fig. 4: (Color online) Partial spin-majority and spin-minority 
Ni 3d DOS (per Ni atom) for spin-polarized La5/3Sr!/ 3 Ni04 in 
the CDFAS phase, as obtained by the LDA+C7 method. 



reference for the total energy gain of specific spin patterns. 
In particular, we investigate diagonal filled ferromagnetic 
stripes (DFFS) as well as A- and C-type diagonal filled an- 
tiferromagnetic stripes (ADFAS, CDFAS), which allows us 
to compare ferromagnetic with antiferromagnetic coupling 
both within the NiO planes (intraplane) and between ad- 
jacent planes (interplane) . 

The electronic states at the Fermi energy are expected 
to grow out of the Ni 3d and 2p orbitals, due to Ni-0 
bonding. This is confirmed by fig. El which shows the Ni 
and O partial densities of states (DOS) resulting from a 
spin-degenerate LBA+U calculation. Contributions of La 
and Sr states in the energy interval of fig. [2] are very small 
and play no role for the further discussion. Here, and in 
the following LBA+U calculations, local interactions are 
assumed for both the Ni 2+ and Ni 3+ sites. For the onsite 
Coulomb interaction we choose a value of U = 8eV and 
for the Hund's coupling a value of Jr = 0.8 eV [28,29]. 
In fig. [2] we identify two groups of electronic bands, the 
first extending from about — 8eV to — 1.6 eV, with respect 
to the Fermi energy. The second, being separated from 
the first group by an energy gap of 0.7 eV, starts at about 
— 1.9eV and extends far beyond the Fermi level. Whereas 
2p contributions concentrate in the low energy group 
of bands, the DOS at higher energy is dominated by the 
Ni 3d states. In particular, a pronounced DOS structure 
is evident at the Fermi energy, which points towards an 
instability against magnetic ordering. 

Introduction of the spin-polarization has serious effects 
on the DOS of the magnetic Ni atoms. In figs. [3] and [4] 
we address the spin-polarized partial Ni 3d DOS, sepa- 
rated into contributions from the four crystallographically 
inequivalent Ni sites, for two possible antiferromagnetic 
spin patterns of La5/3Sri/3Ni04. The magnetic structure 
of the ADFAS phase (fig. [3]) is based on a ferromagnetic ar- 



rangement of domain walls and magnetic domains, an in- 
traplane ferromagnet, and antiferromagnetic coupling be- 
tween adjacent planes. On the contrary, in the case of the 
CDFAS (fig. [4]) the diagonal stripes are ordered antifer- 
romagnetically within the NiO planes, as depicted in fig. 
[U whereas the coupling is ferromagnetic along the per- 
pendicular c-axis. For both antiferromagnetic spin pat- 
terns, the DOS shows a broad structure of occupied Ni 3d 
states, and a finite number of states at the Fermi energy. 
Lacking of an insulating band gap, however, contradicts 
an experimental gap of 0.26 eV [19,20]. In figs. [3] and |4] 
we find broad spin-minority bands extending far beyond 
the Fermi level. For the ADFAS phase the same applies to 
the spin-majority states, while for antiferromagnetic intra- 
plane coupling these states give rise to a much more pro- 
nounced structure with a remarkable weight at the Fermi 
level. In addition, we have tried several other spin pat- 
terns in order to obtain an insulating state, but without 
any success. In each case, the system appeared to be far 
from opening a band gap. Furthermore, we have investi- 
gated whether a local interaction on the oxygen sites has 
to be taken into account. While the total energy in fact 
slightly decreases for reasonable values of U and Jh, the 
DOS shape changes only marginally near the Fermi level. 
In a recent LDA+/7 study on La5/3Sri/3Ni04, Yamamoto 
et al. [30] have found an insulating solution with an energy 
gap of 0.11 eV. However, in contrast to our calculations, 
their data rely on a spin pattern with antiferromagnetic 
ordering along the diagonal stripes. 

We next address the stability of the spin patterns under 
consideration. According to tableEl all magnetic solutions 
are lower in energy than the paramagnetic phase, already 
for U — 0. In particular, the DFFS and ADFAS phases 
turn out to be almost degenerate, both revealing an energy 
gain of 5.8mRyd per Ni atom as compared to the non- 
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DFFS 


ADFAS 


CDFAS 


ADFAS, LBA+U 


CDFAS, LBA+U 


intraplane coupling 


fe 


fe 


af 


fe 


af 


interplane coupling 


fe 


af 


fe 


af 


fe 


magn. moment Nil (/^b) 


0.72 


0.66 


1.14 


1.66 


1.73 


magn. moment Ni2 (/ib) 


0.71 


-0.71 


1.10 


-1.53 


1.63 


magn. moment Ni3 (/is) 


0.77 


0.72 


-1.02 


1.68 


-1.60 


magn. moment Ni4 (/ie) 


0.78 


-0.78 


-0.89 


-1.64 


-1.55 


total energy gain (mRyd) 


5.8 


5.8 


87 


152 


143 



Table 2: Ni magnetic moments and total energy gain (per Ni atom) with respect to non-magnetic calculations. Columns 1 to 3 
refer to pure GGA calculations, whereas columns 4 to 6 give the corresponding LDA+(7 values. The LDA+C/ results are based 
on the parameters U = 8 eV and J=0.8 eV. 



magnetic solution. Because these patterns only differ by 
their spin arrangement along the c-axis, we conclude that 
the magnetic interlayer coupling is very weak. Turning to 
the comparison of the ADFAS and CDFAS configurations, 
we obtain that the latter has a substantially lower energy 
than the former (6.7 mRyd per Ni atom). This energy gain 
can safely be attributed to the intraplane magnetic cou- 
pling, due to the nearly vanishing interplane interaction, 
as mentioned above. Remarkably, the evaluation of the lo- 
cal magnetic moments self-consistently influences the var- 
ious hopping matrix elements and crystal fields involved 
in the calculation. It therefore is unlikely to design a re- 
alistic model without accounting for this self-consistency. 
Nevertheless, a large number of calculations for the two- 
band Hubbard model clearly point at competing phases, 
see [31,32] and the references therein. 

Regarding the magnetic moments, we find comparable 
amplitudes for all Ni sites, independent of the magnetic 
structure. Remarkably, even in the case of the DFFS con- 
figuration the magnetic moments show a strong tendency 
to stripe formation. Their values increase substantially in 
the CDFAS phase, see table resulting in a lower energy. 
Note that the total magnetic moment is finite in one unit 
cell, amounting to 0.35 [Ib per Ni site. Nevertheless, the 
moment is expected to vanish if the supercell would be 
doubled in the c-direction. Since the interlayer coupling is 
almost negligible, the result of such a calculation would be 
closely related to the present data, yet producing a van- 
ishing total magnetization. Our calculation thus captures 
all relevant features of the CDFAS phase. Moreover, a so- 
lution with vanishing intraplane magnetic moment could 
not be obtained, and therefore is likely to be either higher 
in energy or not to exist. 

Concerning the magnetic state of the Ni ions, we find 
a superposition of S z — 1/2 (Ni 3+ -like) with S z = Q and 
S z — 1 states (Ni 2+ -like, one electron in each e g orbital 
with parallel/antiparallel spin), whereas states involving 
doubly occupied orbitals are suppressed [33]. When the 
onsite electron-electron interaction is included (U = 8 eV, 
Jh = 0.8 eV), the local magnetic moments increase signif- 
icantly, to about 1.6/xg in both antiferromagnetic phases. 
This modification arises from a combined effect of U and 



Jh, as suggested by model calculations [33]. However, the 
increase is amplified in the ADFAS phase as compared 
to the CDFAS phase, favoring the former by 9 mRyd per 
Ni atom. Thus, the energy order of the antiferromagnetic 
patterns is inverted due to the onsite interaction, hence no 
longer reflecting the experimental situation. This fact in- 
dicates strong interplay between the magnetic structure, 
the onsite Coulomb interaction, and the energetical se- 
quence of the magnetic phases. We have verified that the 
latter is not sensitive to (i) the exact values of U and Jh , 
in a reasonable range compared to experimental data, and 
(ii) the inclusion of a local interaction on the O sites. Com- 
bining the above with the absence of an insulating ground 
state in our band structure calculation leaves little room 
for an explanation of the deviation from the experimen- 
tal observation: one appealing scenario is provided by the 
fact that the final determination of the crystal structure, 
in particular the positional parameters of all atoms, has 
not been performed so far. The structural estimates used 
for the atomic positions can be inadequate, suggesting to 
perform a structure optimization. Up to now, however, a 
realistic structure optimization could not be realized due 
to a large number of lattice degrees of freedom. As struc- 
tural relaxation is expected to be very sizeable for NiC>6 
octahedra, modification of the Ni-0 bonding could easily 
alter the electronic structure in the vicinity of the Fermi 
energy, such that an insulating gap arises. Of course, this 
might correct the energy sequence of the spin patterns. 



In summary, we have presented a series of spin-resolved 
electronic structure calculations for the striped nickelate 
La5/3Sr!/ 3 Ni04, including onsite Coulomb interaction. In 
contrast to experimental findings, the LDA+U approach 
results in A-type diagonal filled antiferromagnetic stripes, 
while in the pure GGA scheme C-type diagonal filled an- 
tiferromagnetic stripes are favored. Our results indicate 
that accordance with the experimental situation calls for 
further determination of the structural parameters, since 
NiOe octahedra are known to be subject to strong struc- 
tural distortions in many cases. We expect that the same 
applies to all commensurate doping levels. 
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